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ABSTRACT 

A model for element abundance fractionation between the solar chromosphere 
and corona is further developed. The ponderomotive force due to Alfven waves 
propagating through, or reflecting from the chromosphere in solar conditions 
generally accelerates chromospheric ions, but not neutrals, into the corona. This 
gives rise to what has become known as the First Ionization Potential (FIP) 
Effect. We incorporate new physical processes into the model. The chromo- 
spheric ionization balance is improved, and the effect of different approximations 
is discussed. We also treat the parametric generation of slow mode waves by the 
parallel propagating Alfven waves. This is also an effect of the ponderomotive 
force, arising from the periodic variation of the magnetic pressure driving an 
acoustic mode, which adds to the background longitudinal pressure. This can 
have subtle effects on the fractionation, rendering it quasi-mass independent in 
the lower regions of the chromosphere. We also briefly discuss the change in 
the fractionation with Alfven wave frequency, relative to the frequency of the 
overlying coronal loop resonance. 

Subject headings: Sun:abundances - Sun: chromosphere - turbulence - waves 

1. Introduction 

The First Ionization Potential (FIP) effect is the by now well known enhancement in 
abundance by a factor of 3-4 over photospheric values of elements in the solar corona with 
FIP less than about 10 eV. These low FIP elements include Fe, Si, Mg, etc. Elements 
with FIP greater than 10 eV (high-FIP) mostly retain their photospheric composition. This 
was actually first observed in the 1960's (Pottasch 1963), making it nearly as old as the 
problem of coronal heating. It has been taken seriously as a phenomenon since the mid 
1980's (Meyer 1985ab). There are a number of studies of the FIP effect in different regions 
of the solar corona and wind, reviewed in Feldman & Laming (2000) and Laming (2004), 
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and references therein. With the launch of the Extreme Ultraviolet Explorer (EUVE) in 
the 1990 's, it became clear that element abundances in late-type stellar coronae also do 
not always resemble the stellar photospheric composition. The FTP effect is also observed 
in many solar-like late type stars. At higher activity levels and/or later spectral types, a 
so called "i nverse FIP" effect is o bserved, where the low FIP elements are depleted in the 



corona (e.g. IWood fc LinskyI 120101 ). 



A variety of models have been proposed to explain these phenomena. iLaming ( 12004 



20091 ) review many of these different scenarios, and argue that the ponderomotive force is 
the most likely agent of FIP fractionation. This force arises as Alfven waves propagate 
through the chromosphere, and acts on chromospheric ions, but not neutrals. Physically, 
it corresponds to the interaction of waves and plasma through the refractive index of the 
medium. Waves carrying significant energy and momentum refracting or reflecting in a 
plasma must exert a force, in this case on the charged particles that contribute to the 
dielectric tensor, but not on the neutrals. The ponderomotive force in the chromosphere 
may in principle be directed upwards or downwards, giving rise to FIP or so-called "inverse 
FIP" effects respectively (i.e. a coronal enhancement or depletion of low FIP ions). 

The chromosphere-corona interface is generally a barrier to Alfven wave propagation; 
upcoming waves from the chromosphere are usually reflected back down again, and down- 
ward directed waves from the corona reflect back upwards, as illustrated in the right hand 
footpoint (chromosphere B) in Figure 1. Alfven waves with predominantly coronal origin 
generally give rise to the positive (i.e. solar-like) FIP effect, while waves generated by up- 
ward propagating acoustic waves associated with stellar convection may produce inverse FIP 
effect. The association of coronal abundance anomalies with Alfven waves gives us a unique 
and unexpected diagnostic with which to explore the behavior of MHD turbulence in solar 
and stellar upper atmospheres. While we will argue that the coronal Alfven waves themselves 
are actually byproducts of processes that heat solar and stellar coronae, (most likely accom- 
plished by various forms of "nanoflares" ) , an understanding of coronal abundance anomalies 
still becomes far more central to exploring coronal heating than would be the case in prior 
models for the fractionation invoking thermal processes such as diffusion. 

In this paper we seek to develop the ponderomotive force model incorporating the para- 
metric generation of parallel propagating slow mode waves by the Alfven waves themselves, 
together with revisi ons to some of the ato mic data. First, in section 2, we review the impor- 



tant features of the iLamingl (12004 . l2009l ) model. Section 3 outlines the various theoretical 
refinements made in this paper, and section 4 describes new results for fractionations in open 
and closed magnetic field configurations. Section 5 discusses in more detail the effect of the 
parametric generation of slow mode waves and its implications for fractionation in closed 
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coronal loops and in the slow speed solar wind. We also consider limits on the upward flow 
speed through the chromosphere, and make final conclusions. 



The Ponderomotive Model Revisited 



The ponderomotive force stems from second order terms in 6 J x 5B/c and p5v ■ V5v in 
the MHD momentum equation. It can be manipulated (e.g Litwin & Rosner 1998, Laming 
2009) for waves of frequency ua << ^, the ion gyrofrequency, into the time averaged form 
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where 6E± is the perpendicular wave electric field and q is the ion charge. The ponderomotive 
acceleration, F/m, is independent of the ion mass. The Laming (2004, 2009) model comes 
about as a natural extension of existing work on Alfven wave propagation in the solar 
atmosphere with essentially no extra physics required. It is also the model most worked out 
in detail to interpret observations, giving it unique potential for diagnosing wave processes 
in the corona and chromosphere. 

The basic model builds on HoUweg (1984), where upward propagating Alfven waves 
were introduced at one loop footpoint. Here they could refiect back down into the chromo- 
sphere, or be transmitted into the loop, where they propagated back and forth with a small 
probability of leaking back into the chromosphere at each end. With reference to Figure 1, 
we initiate our simulations with one downward propagating wave at the /3 ~ 1.2 layer in 
chromosphere A, and integrate the Alfven wave transport equations (see below) to chromo- 
sphere B to evaluate the st anding wave pattern there. The chr omosphere at each footpoint 
can be based on any of the IVernazza. Avrett. fc Loesen (jl98ll ) models or similar. Here we 
use the Avrett Sz Loeser (j2008l ) update of VALC. The ionization balance of the minor ions 
is computed at each height in the chromosphere using the model temperature and electron 
density, and a coronal UV-X-ray spectrum appropriately absorbed in t he intervening chro- 
mosph eric layers. We have tried a number of different spectra bas ed on IVernazza &: Reeves 



(119781 ) ■ or model fiare spectra computed using CH IANTI (see e.g . iHuba et al. 



2000 Bastille Day fiare). The atomic data are as in iLamingl (120041 ) and iLamingl ( l2009l ). with 



2005L for the 



estimates for the charge exchange ionization for Si, Fe, and other low FIP elements added 
(see subsection 3.3). The chromospheric magnetic field is taken to be a 2D force free field 
from Athay (1981) and designed to match chromospheric magnetic fields in Gary (2001) 
and Campos & Mendes (1995), which represents the expansion of the field from the high 
P photosphere where the field is concentrated into small network segments, into the low /3 
chromosphere where the field expands to fill much more of the volume. 
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We model the Alfven waves in a non-WKB approximation. The procedure follows that 
described in detail by Cranmer & van Ballegooijen (2005), but applied to closed rather than 
open magnetic field structures. The curvature of the loop is neglected. For Alfven waves, 
where the energy fiux is necessarily directed along the field direction, this is unlikely to have 
a significant effect. Some extra damping may result as the wave develops a component of its 
wavevector perpendicular to the field, but we neglect this and all other damping mechanisms 
in this work. The transport equations are 

dZ± ^ , , TA \ I , TA \ ^± , % 



^ + (2) 

where z± = Sv± ± SB± / ^/^np are the Elsasser variables for Alfven waves propagating against 
or along the magnetic field respectively, and are valid for torsional or planar Alfven waves. 
The Alfven speed is Va, the upward fiow speed is u and the density is p. The signed scale 
heights are Hd = p/ (dp/dr) and Ha = VaI idVAldr). In the solar chromosphere and corona 
u «Va, and we put m = 0. The calculation of Va uses the total (ionized and neutral) gas 
density, since the wave frequencies of interest here are well below the charge exchange rate, 
and neutrals are well coupled to the wave. Charge changing collisions involving electrons 
(impact ionization, and radiative and dielectronic recombination) are generally slower than 
charge exchange in chromospheric conditions. Hence in considering the wave propagation, 
ion-neutral friction is neglected, though it is included in the evaluation of the fractionation. 

The fractionations are calculated by postprocessing the non-WKB wave. This is valid 
because the fractionation has a negligible effect on the wave propagation. The degree of 
fractionation is given by the formula 

fractionation = exp ^2 j CsO'^eff/^s,i/'Vsdz^ (3) 

(see equation 9, Laming 2009, equation 12, Laming 2004, which follow Schwadron et al. 
1999), where is the ionization fraction of element s, a is the ponderomotive acceleration, 
I \isys,n + (1 — Cs) Vs,i\ where Vs^i and Vg^n are the collision rates of ions and 
neutrals, respectively, of element s with the ambient gas. Also v1 = k T/nis + vf,-t^,rb +. 



Vfyrb^ where Vfj,turb is the amplitude of microturbulence coming from the lAvrett fc Loeser 



(120081 ) chromospheric model, and Vturb, discussed further in section 3.2, is the amplitude of 
longitudinal waves induced by the Alfven waves themselves. The limits of integration, zi and 
Zu are the lower and upper limits over which the ponderomotive force acts. We take zi to be 
where /3 — 1, and z„ is in the transition region at an altitude where all elements are ionized. 
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3. New Model Features 

3.1. Slow Mode Waves 

We have introduced an extra longitudinal pressure associated with the Alfven waves 
proportional to v^^^.^,, which has the effect of causing some saturation of the FTP fractionation. 
Here we give the physical motivation for this extra term, which arises from the generation of 
slow mode waves. Physically, the periodic variation in magnetic pressure of the Alfven wave 
drives longitudinal compressional waves. These generated acoustic waves can act back on 
the Alfven driver, as the compressional wave introduces a periodic variation in the Alfven 
speed, which generates new Alfven waves. We illustrate the generation of slow mode or 
acoustic waves by the ponderomotive force associated with plane Alfven waves with a simple 
ID calculation. The linearized momentum equation keeping terms to all orders in perturbed 
quantities is (all symbols have their usual meanings), 

, - , dSvz , „ , „ dbvz , c- X 9 bB"^ d5P . , 

(P + ^ + (p + = (p + ajg^f^ - ^ - s*/'. (4) 

where 

for Eulerian displacement ^ oc expi {ust + kgz), where Lp — P/ (dP/dz) and Lp — p/ (dp/dz) 
(signed) pressure and density scale heights respectively. The first term on the right hand side 
of equation 4 represents the instantaneous ponderomotive force. In cases where the WKB 
approximation applies, SB oc p^^"^, and this expression is equivalent to the more usual form 
—d{6B'^/8'7i) /dz. Substituting for 5p and 5P from equations 5, keeping terms as high as 
second order in small quantities, equation 3 becomes 

- 'Tf'^ + ^-"^ + ^^'^ -'t;-Wp~'Wp~ -Tj^'^-'^^'d-zs^p^'- 

This is considerably simplified in isothermal conditions, 7 = 1, Lp — Lp — —c^/g, so that 

UJg 2 / 2 t22 ^ SB^ 

- i—ov^ + l-Ug + kgCg + iksg) dv^ - ^cu^— - — = 0. (7) 
Lp ^ ' oz Snp 



We put ^kg — —g/2cl, and y {^ks) + g'^/4:C'^ = 23^^/c^/r^, u — 2ua/ti, where /c^ and uja are 
the wavevector and angular frequency of the Alfven wave with n — 1,2,3, .. (anticipating 
the result below). We find 
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with solution 



5v, 



—I 



— LpbJs 
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(9) 



where we have put ^ {5B'^/8Tcp) = {6B'^/ATcp) /ALp = 5v\/ALp using the WKB resuh for 
Alfven waves in a density gradient, and assuming Ij Lp » 3ft A;^. When ~ V| or Lp — i- 0, 
5vz — —idvA/'^- In the opposite hmit 5vz — —i5v\/ ALpUJs (1 — c^/V|)- In these two cases 
ojs = (jJa or ojs = 2uja respectively. In fact acoustic waves can be ge nerated with 0;.^ = 2u)A/n , 
with higher n becoming mo r e inte nse as the nonhnearity increases (iLandau &: Lifshitall976l ). 
Vasheghani Farahani et al.l (I2OIII ) treat the case of slow mode wave generation by a torsional 
Alfven wave. This is subtly different to the case of a plane Alfven wave considered here, and 
the FIP fractionation resulting from such a model will be investigated in a future paper. 



Anticipating applications to possibly nonlinear Alfven and slow mode wave amplitudes, 
we revisit the analysis above retaining more higher order terms. From 6p = —pC, {ikg + 1/L) 
we derive 



which when substituted into the hnearized continuity equation. 
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with the time derivatives d5vz/dt = iuJsSvz and d5p/dt = iug^p, gives 
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p5vz 



1 + 



2k Jv. 



+ 



SVz 

iuL, 



:i2) 



Writing 6P = 7P {Sp/ p + 6vz/iujLp) — P6vz/iujLp we similarly derive 
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We now eliminate Sp and 6P in favor of 6vz in equation 3 to derive a quartic equation in 6vz 
for the driven slow mode wave with angular frequency Us and wavevector ks', 
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0. 



(14) 



To lowest order, the terms in 6vz and the constant are the same as in equation 6. The 
quadratic and higher terms are changed, because of the difference between equations 12 and 
13, and those in equation 5. Inserting even accurate spatial derivatives in place of those in 
equations 12 an d 13 would generat e yet m ore higher order terms, extending n in principle 
without limit. iLandau fc Lifshitzl (jl976[ ) give a similar conclusion in their treatment of 
parametric resonance. 

We solve equation 14 numerically, with the same kg (real and imaginary parts) as above. 
We select the solution with the lowest absolute magnitude as the physically correct solution 
for our problem, this being the solution that goes to zero as 6va — ^ 0. This is always close 
to the solution obtained discarding all terms of order higher than Sv"^ in equation 14, and 
usually close to the case when the term in 6v't is also ii e glecte d. This describes turbulence 
for parallel propagating waves. IZaqarashvili fc Roberts! ( 120061 ) give a detailed treatment of 
the interaction between weak Alfven and sound waves in a homogeneous medium, where 
acoustic and Alfven speeds are equal. The stronger generation of acoustic waves by the pon 



deromotive force in a density gradient is demonstrated by the simulations of iDel Zanna et al. 



teOOSi ). where slow mode waves are seen propagating up from loop footpoints with properties 
consistent with the solutions of equations 9 or 14. 

Closer to the layer where = V|, the magnetic field becomes more curved, giving rise 
to higher perpendicular components of Alfven wave wavevectors, and potentially stronger 
tu rbulent c a scade and/or mode conversion. In this case we expect stronger slow mode waves. 
In iLamingl (j2009[ ) we assumed 6vs = Sva- However considerations of mode conversion for 
initially upward propag ating acoustic waves sugge st that higher slow mode intensities than 
this should be present. iKhomenko fc Callyl ( 1201 ll ) report that at conditions for maximum 
conversion of a high /3 fast mode wave to a low /3 Alfven wave, the low /3 slow mode wave 
has 2-3 times more flux. In this paper, we make the approximation 



6vi = 6vi + 65vici/V, 



r2 

'Ay 



(15) 



where 6vz represents the solu tion of equation 14, an d th e factor 6 in t he las t term is motivated 
by calculations illustrated in lCranmer et al.l (120071 ) and iKhomenkd (|2010[ ). Slow mode waves 
governed by equation 15 have the effect of suppressing fractionation in the low chromosphere 
close to where the plasma /3 ~ 1.2. Studies of mode conversion between acoustic and Alfven 
waves generally show the upward moving acoustic wave beginning to convert to Alfven waves 
at the P ^ 1.2 layer, and mode con version continuing over a range of altitudes of order lOO's 
of km (e.g lCally fc GoossensI 120081 ). The explicit incorporation of such effects is well beyond 
the scope of the work here, and the prescription in equation 15 should be sufficient to avoid 
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the occurrence of unphysical fractionations. Even so, it remains a feature of this work 
requiring further investigation in future papers. 



3.2. Normalization Relative to Oxygen 



In previous papers iLamingl (120041 . |2009| ) we have discussed FIP fractionations relative 
to H. However the derivation of 3 has assumed fractionated elements are minor species, with 
the fractionation having no back reaction on the flow due to the neglect of inertial terms. 
Thus it is more appropriate to present and describe element fractionations with respect to 
another minor element. We choose O, which is a common choice for observers also. 



3.3. Charge Exchange Ionization 



Charge exchange rates have been previously taken from the compilation of lKingdon fc Ferland 



( 119961 ). These have been supplemented more recently with charge exchange ionization rates 
for Si and Fe, colliding with protons. We implement an estimate of the charge exchange 
ionization rate for all ions with lower FIP than H as follows. We estimate the radius R at 
which the sum of the binding energy of the electron in the target neutral and the polarization 
potential energy of the target neutral in the electric field of the incoming proton and equal to 
the equivalent sum for the resulting neutral H atom in the electric field of the newly formed 
ion, (known as the radius of the potential crossing) from 



V 



[as - aH) 



-Ih + Is 



(16) 



where V represents the difference in potential energy of the proton in terms of the polar- 
izability of the target atom and resulting ion in terms of an, the polarizability of the 
resulting hydrogen atom. Ih and Is are the ionization potentials of hydrogen and the target 
atom, s, respectively. Polarizabilities and ionization potential here are in atomic units. The 
cross section is then acxi = t^R? I'2. = 7ry^(a7^-"tt/;T7T^^r--Xy/2, assuming the maximum 
probability for a reaction is 1/2 . Thi s estimate is a slightly different form of the Langevin 
formula given bvlFerland et al.l (119971) . This approximation gives good agreement with the 
calculations of Allan et al. ( jl988 ) for charge exchange ionization of Mg. These authors com- 
ments that similar rates should exist for all other elements with ionization potentials below 
that of H, and we apply it to all of these elements. The charge exchange implemented is 
just the thermal process. No account is taken of any possible effects of the waves on the 
chromospheric ionization balance. 
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3.4. Comparison of Different Approximations 



Carlsson fc Steiru (20021) have argued that the concept of an "average" chromosphere 
pursued by lAvrett &: Loeserl (120081 ) and its antecedents is invahd, due to the extreme dy- 
namics associated with chromospheric shock waves, arguing that "the mean value of a 
dynamic property is n ot the same as that property evaluated for the mean atmosphere". 
Avrett fc Loeserl (120081 ) derive mean values for plasm a properties based on ob servations, not 
on calculated mean atmospheres. More importantly, ICarlsson fc SteinI (120021 ) show that the 
ionization fraction for H varies very little about its mean, due to the length of ionization and 
recombination times compared to the fre quencies of shocks, and that t heir a verage electron 
density agrees very well with that in the IVernazza. Avrett. fc Loeserl ( 198ll ) model C. It is 
easy to see that other high FIP elements should behave similarly, and that the ionization 
balance we calculate (the overwhelmingly most important chromospheric property to us) 
should not be greatly in error, if at all. The inclusion of charge exchange ionization (sec- 
tion 3.3) increases the ionization level for all other elements as the ionization of hydrogen is 
increased above its thermal equilibrium l e vel. T he Ca"*" to Ca^+ ionization balance is con- 
sidere d bv IWedemever-Bohm fc CarlssonI (120111 ). This is more variable than that for H to 



H"*" in lCarlsson fc SteinI ( 120021 ) . but is less of a concern to us, since the ponderomotive force 



experienced by an ion is independent of ion charge, so long as uja « ^■ 

We give sample calculations that go some way to quantifying the effect of these issues. 
Figure 2 shows the coronal section of the 100,000 km long loop with magnetic field B = 20G. 
The density at the loop apex is 5 x 10® cm~^. This gives a resonant angular frequency 
of 0.07 rad s~^. A wave of this frequency propagates on the loop, which is thus half a 
wavelength long. The top panel gives the Elsasser variables (real and imaginary parts), 
the middle panel gives the wave energy fluxes and their difference, and the third panel 
gives the ponderomotive acceleration. Figure 3 shows the chromospheric section of the 
loop where the fractionations are evaluated. This has the same three panels as for Figure 
2, where the third panel also gives the slow mode wave amplitude, with a fourth panel 
(bottom right) that gives FIP fractionation and the ionization fraction of elements Fe, Mg, 
S, and He relative to O. In Table 1 fractionations for He, C, N, Ne, Mg, Si, S, Ar, and 
Fe relative to O are displayed, calculated according to our basic model described above, 
labeled "baseline". In the succeeding columns, we give FIP fractionation s calculated with 
differe nt modiflcations to the model. In the flrst case the density given by lAvrett fc Loeser 
J2008h is modifled so that the model density is consistent with the H ionization fraction and 
the assumption of photoionization-recombination equilibrium. The second variation gives 
fractionations calculated assuming the ionization fractions to be given by the Saha equation 
at the temperature and density in the chromospheric model. The flrst variation reduces the 
degree of ionization in the chromosphere, while the second one increases it. This has the 
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opposite effect on the fractionations, since tliese are given relative to O, and tlie increase 
or decrease of tfie ionization of O lias a bigger effect on its fractionation than is the case 
with the other elements. As can be seen, the basic phenomenon of the FIP effect remains 
unaffected by the choice of approximation, but some of the details, e.g. the Mg/Fe ratio are 
subtly different. We return to discuss this in more detail below, in subsection 4.2. 

The final column in Table 1 gives the FIP fractionation calculated with the assumptionof 
slow m ode wave amplitude 5vz = Sva, the amplitude of the Alfven wave, as taken in iLaming 
( 120091 ). instead of implementing the solution of equation 14. The overall fractionation is 
reduced for the same Alfven wave amplitude, due to the inc rease in longitudinal pressure 
with the higher slow mode wave amplitude. In iLamingl (j2009[ ).we found that the FIP Effect 
saturated in this case at values around 3-4 for arbitrarily high Alfven wave amplitude. In 
section 4, we will consider the behavior of the FIP Effect with resonant and nonresonant 



waves. 



Results 



4.1. Coronal Hole 



A coronal hole is modeled int he same fashion. The chromosphere is identical to the 
case above, and the density evolves smoothly o ff- limb, declinin g to about 2 x 10^ cm~'^ at 
an altitude of 5 x 10^ km (as in e . g. Fig ure 5 of ICranmerl |2009| ) . The magnetic field follows 
the model of 
calculated by 



Banaskiewicz et al.l (I19981) . We take the coronal hole Alfven wave spectrum 



Cranmer et al.l ( 120071 ) as our starting point. It is illustrated in their Figure 3 



for a position in the solar transition region. It is represented by the five wave frequencies 
and amplitudes at the starting position of the non-WKB integration of the wave transport 
equ ations, in this case at a n altitude of 500,000 km. Parameters are chosen to match Figure 
9 in ICranmer et al.l (120071 ). and are given in Table 2 as the spectrum labeled "vO". Figure 
4 shows the wave amplitudes in the coronal hole section of the calculation, with the same 
three panels as for Figure 2, but illustrated up to an altitude of 500,000 km. Figure 5 shows 
the chromospheric response with the same four panels as in Figure 3. The ponderomotive 
acceleration (Figure 5, top right) has a similar "spike" at the top of the chromosphere as in 
Figure 3 for the closed loop, but is about an order of magnitude weaker. Lower down, the 
Alfven wave amplitudes, and the corresponding slow mode wave amplitudes are larger than 
before. The net effect of smaller ponderomotive acceleration and higher slow mode wave 
amplitude is to reduce the FIP effect from the case i n Figu re 3. The resulting fractionation 
here is very similar to that found by I Cranmer et al.l ( l2007l ). In models vl-2 in Table 2, we 
increase the amplitude of the highest and lowest frequency wave in the spectrum respectively. 
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Increasing the highest frequency wave amphtude has the effect of increasing the fractionation 
predicted, while increasing the low frequency wave has very little effect. The results from 
all three models are summarized in Table 3, and compared with observations for various 
open field regions. Reasonable agreement for the fractionation of low FIP elements is seen. 
A small depletion of He is seen, but not as large as observed. Our purpose here has not 
been to provide a definitive calculation of the FIP effect in a coronal hole, but merely to 
illustrate that the observed difference between the FIP effect in closed and open field lines 
arises naturally in this model. 



4.2. Closed Loop 



We now explore fractionation in a closed coronal loop, illustrating the difference that a 
coronal resonant frequency can make to the fractionation. In Table 4 we give the fraction- 
ations computed for a closed loop with length 100,000 km, and a magnetic field of 20 G. A 
coronal wave at the loop resonant angular frequency of 0.07 rad s~^ is modeled. The wave 
transport equations (equation 2) are integrated from the /3 = 1 layer in one chromosphere, 
where and initially downgoing wave amplitude is specified, through the loop to the opposite 
chromospheric footpoint, where the FIP fractionations are evaluated. The initial wave am- 
plitudes are 0.4, 0.5, and 0.6 km s~^, giving coronal wave peak amplitudes of approximately 
55, 70, and 82 km s~^ respectively. 

These cases illustrate the variation of the FIP effect with wave amp litude. The six 
columns on the right hand side give v arious observations of FIP fractionation. IZurbuchen et al 
( I2OO2I ) and Ivon Steiger et al.l (I2OOOI ) both give fractionations mea s ured i n situ in the so- 
lar wind over relatively long periods of time. iGiammanco et al.l (120071 120081 ) give frac- 
tionations also measured in the solar wind, but ove r time periods select e d such that the 
wind s pee d was close to 380, 390 or 4 00 km s~^ 
feoiOal Jbh: iMcKenzie fc Feldman 



Phillips et al.l ( 120031 ): ISylwester et al. 



■ ■ _ 19941) give abundance ratios measured spectroscopically 

in solar flares, and Bryans et al. ( 20091 ) observe quiet solar corona above the western limb. 



also spectroscopically. While the agreement between the FIP fractionations calculated for 
the initial Alfven wave amplitude corresponding to 0.5 km s'^ is generally good, there are 
some important discrepancies. The ratio C/0 is typically observed in the solar wind to be 
higher than calculated, as is S/0 for some of the observations. While Fe/0 and Mg/0 are 
reasonably w-ell repr oduced, the direct ratio between Fe and Mg is not, except in the case of 
Bryans et al.l (120091 ). The last case, with initial wave amplitude 0.6 km s~^, is designed to 



match these observations, and does quite well. Only K is seriously discrepant, with S also 
somewhat underestimated. In the solar wind and flare observations, Fe and Mg are often 
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fractionated by the same amount. 

Table 5 gives FIP fractionations for varying Alfven wave frequency, with the amphtudes 
chosen to keep the Fe/0 abundance ratio close to 4 as is often observed. Here, we have also 
included a wave field of chromospheric origin designed to have the same properties as that in 
the vO model for the open field case. Downgoing wave amplitudes are specified at the /3 = 1 
layer in one chromosphere, and then the wave transport equations are integrated back to the 
opposite footpoint where FIP fractionations are evaluated. The column corresponding to 
the resonant frequency of 0.07 rad s~^ is the same as in Table 4, but for the new model. The 
chromospheric waves can be seen to have rather little effect, with the biggest changes being 
seen in the increased fractionations of C/0 and S/0. Moving away from this resonance, 
either to lower or higher frequency, the C/0 and S/0 ratios increase to better agree with 
observations. At higher frequency so too does the Mg/ O ratio, so that Mg and Fe fractionate 
more closely to the same degree; quasi-mass independent fractionation is achieved. Going to 
yet higher frequency, all high FIP elements remain unchanged, while all low FIP elements are 
en hanced by a f a ctor 3 -4. The depletions of He and Ne are lost. The quiet sun observations 



of iBryans et al.l ( 120091 ) are best matched at the frequency of 0.075 rad s~^ and with the 
exception of K, the consistency between their measurements and the model is excellent at 
an initial wave amplitude of 0.7 km s~^. 

Figures 6 and 7 illustrate the ponderomotive acceleration and FIP fractionation within 
the chromosphere for the cases of frequencies 0.06, 0.07, 0.085 and 0.105 rad s~^. At 0.06 
rad s~^, the ponderomotive acceleration is positive from about 800 km up, and has a "spike" 
at about 2150 km, with maximum close to 10^ cm s~^. Fractionation of Fe, Mg, and S is 
similar low down, but the fractionation hierarchy Fe > Mg > S is established in the range 
1500-2000 km. For 0.07 rad s~^, the ponderomotive acceleration has a similar, but slightly 
larger maximum at about 2150 km. In response to this Fe, Mg and S undergo a similar 
and small inverse FIP fractionation up to 1500 km, giving way to positive FIP higher up. 
The fractionation pattern Fe > Mg > S is even stronger here than for 0.06 rad s~^, and is 
mainly occurring at the "spike" in the ponderomotive acceleration. In the last two cases, 
the ponderomotive force is stronger lower down in the chromospheric, and the "spike" at 
2150 km becomes less pronounced. The slow mode wave amplitude is also stronger lower 
down. At 0.105 rad s~^, all fractionation occurs by 1600 km, and the local maximum in the 
ponderomotive acceleration at 2150 km has no effect. 

These fractionation patterns have simple qualitative explanations. First we display in 
Figure 8, the 0.07 rad s~^ case again to illustrate the relation of the fractionation to important 
features of the chromosphere. The left panels give the ponderomotive acceleration (bottom) 
and the density and temperature structure of the chromosphere (top). The "spike" in the 
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ponderomotive acceleration can be seen to stem from the steep density gradient beginning 
at an altitude of about 2100 km. The solid and dashed lines in the bottom plot show 
the ponderomotive acceleration with and without the energy loss to slow mode waves. In 
the regions where significant fractionation occurs, the slow mode wave do not affect the 
ponderomotive acceleration very much, and their main effect on the fractionation is through 
the additional longitudinal pressure they provide. The panels on the right hand side show the 
same fractionations as before (bottom), and on the top an expanded view of the ionization 
fraction of the elements C, S, Mg, and Fe. 

With reference to equation 3, in regions where H is predominantly neutral, (below 
about 1500 km altitude) Us^i ~ Ug^n and similar fractionation results for elements where 
is reasonably close to unity. Where H is predominantly ionized, Ug^i » i's,m and small 
departures in from unity can make a big difference to the fractionation. This is the 
reason why S fractionates similarly to Fe and Mg in the low chromosphere, but markedly 
less so in higher regions. This is also the reason why Mg fractionates less than Fe higher 
up. At an altitude of 2000 km, the charge state fractions of Fe, Mg, and S are 0.9995, 
0.9981, and 0.9942 respectively (see Figure 8). Even though these are close to unity, the 
differences from unity result in different fractionations where the H ionization fraction (which 
closely follows that of O) is about 0.6. Lower down, where H is mostly neutral, the different 
ionization fractions matter much less in the fractionation. Recalling the results calculated 
using the Saha approximation for the ionization fractions, we can now see why Mg and Fe 
fractionate much more similarly in this case. The ionization fractions at 2000 km altitude 
are now 0.999974, 0.999995, and 0.9983 respectively, for Fe, Mg, and S. These are much 
closer to unity than before, so Fe and Mg now fractionate to a more similar degree. This 
is to be expected, since the assumption of LTE in the Saha equation suppresses radiative 
recombination rates, since the photons so produced cannot escape, and so the Saha ionization 
fractions will be higher than a more realistic calculation would predict. 

However in Table 1 using Saha equilibrium, Fe and Mg do not behave exactly identically. 
Recalling equation 3, 

fractionation = exp I 2 [ — ^ , --^ ; — :rdz\ (17) 

and remembering that Vturb is the amplitude of slow mode waves generated by the Alfvcn 
waves themselves, we can see that when Vturb and f ^jturb dominate over the ion thermal speed 
(usually Vturb > v^turb), the mass dependence disappears from this part of the equation, 
and will only reside, if at all, in the collision frequencies. In fractionation occurring high in 
the chromosphere associated with the "spike" , where the plasma temperature is increasing 
rapidly up to coronal values, this condition may not be met and mass dependent fractiona- 
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tion can occur. In fractionation occurring lower down near the chromospheric temperature 
minimum, for example in the 0.085 rad case, this condition is met, and quasi-mass 
independent fractionation results. 



5. Discussion 

5.1. The Effect of Slow Mode Waves 

The parametric generation of slow mode wave is a crucial part of the fractionation 
process by ponderomotive forces. One important effect has been to render the fractionation 
quasi mass independent as is often observed. This is demonstrated most clearly in Figure 7a, 
corresponding to a 0.085 rad s~^ Alfven wave, the case which also has the highest slow mode 
wave amplitude, staying close to 10 km s^^ for large regions of the chromosphere where 
thermal speeds are ~ 1 km s~^. Fractionation occurring at the top of the chromosphere 
in the location of the "spike" in the ponderomotive acceleration often retains some mass 
dependence, because the plasma temperature is increasing rapidly here while the slow mode 
wave amplitude is usually small. 



In the case that 5vs ~ 5va as in iLamina (120091 ). the increased slow mode wave ampli- 



tude relative to this work suppresses all FIP fractionation except when the wave frequency 
coincides precisely with the loop resonance, and then all fractionation occurs at the loop top 
and hence is mass dependent. For reasons we discuss below, this coincidence is probably not 
realized ubiquitously in the solar corona. Moreover the assumption of isotropic turbulence 
probably requires a well develope d turbulent cascade, w hich is unhkely to develop with purely 



parallel propagating waves (e.g. iLuo fc Melrosell2006l ). Lower down in the chromosphere as 



the plasma beta approaches unity, the magnetic field becomes more concentrated in network 
segments. The increased curvature of field lines will lead initially parallel propagating waves 
to develop perpendicular components to their wavevectors, and hence turbulent cascade or 
mode conversion become more likely. Our equation 15 above is an attempt to capture this 
behavior, and obviously needs to be revisited with greater rigor. 



5.2. The Alfven Wave Frequency 

Table 5 displays FIP fractionations for a range of frequencies close to the fundamental 
of a loop with length 100,000 km and magnetic field 20 G, with wave amplitudes chosen such 
that the fractionation of Fe/0 is close to the usually observed value of 4. We commented 
above how the fractionation details of other elements vary slightly as the coronal wave moves 
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from being in coincidence with the loop resonance, to a position well off resonance. This 
arises because resonant waves reflect from the top of the chromosphere, and this is then the 
sole location of FIP fractionation, but nonresonant waves penetrate further down, allowing 
FTP fractionation to occur over a greater range of altitudes in the chromosphere. 

When FIP fractionation is concentrated at the top of the chromosphere, the different 
ionization structures of the various high FIP elements becomes important, and fractiona- 
tion occurs among them. Most significantly. He becomes depleted relative to O, with this 
depletion being strongest for a frequency 0.075 rad s~^, just higher than the resonance, at 
a value of 0.60. This gives ar i abundance close to the He abundan ce frequently observed 
in the slow speed solar wind (lAellig et al.ll200ll : iKasper et al.l 120071 ). In this frequency re- 
gion too, C and possibly S also have minima in their fractionations. These elements have 
ionization potentials of 11.26 and 10.36 eV respectively (on the boundary between low FIP 
and high FIP elements). Although they are highly ionized throughout the chromosphere, as 
described above, they have sufficient neutral component that they fractionate well when H 
is predominantly neutral, but not when H is ionized. When fractionation is restricted to the 
top of the chromosphere where H is ionized, they behave more like high FIP elements. This 



is commonly seen in spectroscopic observations of S (e.g . iLaming et al.lll995l : iFeldman et al. 
2OO9I : IWiding fc Feldmanll2008l : ISrooks fc Warrenlboilh . 



Off resonance, when FIP fractionation can occur over a more extended range of heights, 
including those where H is mainly neutral, C and S might be expected to behave more 
like low FIP elements. Suc h behavior is more app arent in the solar wind observations of 
Zurbuchen et al.l (120021 ) and lvon Steiger et al.l (120001 ). Here, the FIP bias is variable, so that 
the time average over an extended period gives Fe/0 ~ 2 instead of ~ 4 as modeled. Even 
so, S/0 has a similar value to Fe/0. These observations are best matched in Table 5 for a 
frequency 0.085 rad s~^. 

We have previously suggested that Alfven waves of coronal origin probably derive from 
coronal heating mechanisms such as nanofiares or Alfven resonance. The coronal Alfven 
amplitudes required above (~ 50 — 100 km s^^) are larger than nonthermal mass motions 
observed through spectral linewidths by factors 2-3. This suggests that the Alfven wave 
must be confined to a small fraction of the loop cross-sectional area, which would also be a 
natural consequence of nanofiare or Alfven resonance heating. 

In as far as the heating can be considered a small perturbation to the coronal loop, 
the waves so generated should be eigenfunctions of the loop, with frequencies constrained to 
coincide with the loop resonance(s). The fact that many observations are better matched 
by Alfven wave frequencies slightly higher than the resonance possibly suggests a dynamic 
system. The loop releases waves at its resonance as part of the heating process. The heat 
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conducts dow n to the chromosphere, and heat ed plasma there evaporates back upwards into 



the loop (e.g. iPatsourakos fc KlimchukI |2006| ) . thus increasing its density and reducing the 



coronal Alfven speed, and hence also reducing the loop resonance frequency. The waves 
produced at the original resonance continue to propagate until damped. So we might nat- 
urally expect a mismatch between Alfven wave frequency and the loop resonance, and also 
expect this mismatch to become larger in more strongly heated region of the corona, e.g. 
active region and flares, as opposed to quiet sun. Of course this discussion presupposes that 
the heating is a weak perturbation to the coronal loop, and so its eigenfunctions are well 
defined, and can be excited. In flares and CMEs, this might no longer be true, and the 
heating mechanism itself will determine which waves are produced, irrespective of the loop 
boundary conditions. 

Such ideas will be investigated in greater detail in a separate paper. For the time 
being, we restrict ourselves to some simple predictions. The coronal helium abundance 
shoul d increase with increasing solar activit y, as it appears to do both in solar wind observa- 



t ions (lAellig et al.l 



2001 



Kasper et al.ll2007l) and in spectr oscopic measurement s of the quiet 



sun ( Laming fc Feldman boOl . 2003 ) compared to flares ( Feldman et al. 2p05[). The S and 
C abundance should also vary. In the solar wind (e.g. von Steiger et al. 200ol) it appears 



Laming et al 



1995 



to vary as a low FI P ion, whereas in spectroscopic observations, (e.g. 
Feldman et al.ll2009l . of quiet and active regions) sulfur is observed to behave as a high FIP 
element. 



5.3. The Upward Flovi^ Speed 

We have suggested conduction driven chromospheric evaporation as the source of the 
plasma upflow that populates the corona with fractionated gas. Here we estimate the flow 
speed in the chromosphere, and show that it is consistent with limits set by the operation 
of the ponderomotive force producing the FIP fractionation. With dp/dt = we write 

d , , dp agz dT 

where the density p oc expl—pgz/k^T) is a gravitationally stratifled solution. The mean 
molecular mass is p, and ks is Boltzmann's constant. Integrating between zi and Zu 

p[Zu)v{zu) - p{zi)v{zi) = - j -j^^P'^^- (19) 

We choose the upper limit Zu to be where v {zu) = in the corona, and so 

r-zp(z)d\nT, z ( z\d\nT, , , 
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where Lp = k-BT/fig is the density scale height. This integral will be dominated by the 
integrand near z = zi, so 

where we have put 2.hnk^dT/dt = -V ■ Q with the heat flux Q = -l^'^T^/'^VT . The 
temperature gradient has been replaced by AT/Lt-, where AT may be the coronal peak 
temperature, and Lt the loop half-length. Taking the chromospheric temperature T ~ 10^, 

, , 10^3 / T,/5 X 10^ K V 1 



n \Lt/5 X 10^ cm 

which suggests a velocity of 10^ cm s~^ at a density of lO^'' cm~^. This may well be an 
underes ti mate due to our approximation for the temperature gradient, but as discussed in 



Laming (12004 ). is sufficiently high that gravitational settling should not occur. If v (zi) 



approaches ~ 1 km s ^, as might happen in flares, some further discussion is required. 

The derivation of equation 17 neglected inertial terms in the momentum equations for 
ions and neutrals. Reinstating these, in the limit that Usi — Usn << Ug ^ u the fractionation 
becomes 

fractionation = exp I 2 / -j- — — ^dz 




+ (1 - a i^s,i] [kT/m, + v%^^, + vl^, + u 



flow] 



(23) 

When the magnitude of uj^^^ approaches those of the other terms in the second square 
bracket in the denominator of the integrand, some reduction in the FIP effect will result. This 
is most likely to have an impact on fractionation occurring at the top of the chromosphere, 
taking pujiow ~ constant through the chromosphere, and will possibly reduce the amount 
of mass dependent fractionation occurring there. This will happen for relatively large flow 
speeds Ufiow ~ 1 km s~^ or larger at the top of the chromosphere, still significantly lower 
than the Alfven speed in this region. 



5.4. Conclusions 

We have further developed the model of element fractionation to give rise to the FIP 
effect by the ponderomotive force, paying careful attention to the generation of slow mode 
waves by the primary Alfven oscillations. When considering a realistic wave spectrum with 
both chromospheric and and coronal contributions, the extra longitudinal pressure due to 
the slow mode waves is crucial in producing the correct fractionation. With this extra 
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ingredient, together with the improvements to the ionization balance and the normahzation 
of the fractionation, a rather comprehensive description of the coronal fractionation has 
emerged. In seeking to understand the FIP effect as usually described, we have also found 
an explanation for the depletion of He in the solar wind, and also possibly its variation. The 
Ne abundance also appears to vary in a similar manner, but to a lesser degree. It is also 
more sensitive to assumptions about the ionization balance, but further investigation of this 
is expected to resolve current controversy surrounding the solar photospheric abundance of 
this element. 

The theory now appears to be developed to the point where variations in the FIP 
fractionation from place to place in the solar corona or wind may now be interpreted in 
terms of their physical origins. The element abundances in the solar corona may therefore 
be considered as diagnostics of the behavior of MHD turbulence, and also thereby of the 
mechanisms that heat the solar corona. These ideas will be further developed in subsequent 
papers. 

This work was supported by NASA Contracts NNH10A055I, NNH11AQ23I, and by 
basic research funds of the Office of Naval Research. I am also grateful to Cara Rakowski 
for a critical reading of an earlier draft of this paper. 
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Table 1: FIP Fractionations in Different Approximations 



r3,tio 


HQ cpl irip'^ 


rlpT-icifv morl ^ 


S^^^ h PI inn i f ion ^ 


uUz — u V A 


He/0 


0.67 


0.83 


0.73 


0.85 


C/0 


0.99 


1.26 


1.17 


1.03 


N/0 


0.82 


1.02 


0.95 


0.93 


Ne/0 


0.74 


0.93 


0.90 


0.89 


Mg/0 


1.98 


2.52 


2.33 


1.43 


Si/0 


1.89 


2.37 


2.41 


1.41 


S/0 


1.40 


1.75 


1.47 


1.23 


Ar/0 


0.92 


1.16 


1.03 


0.97 


Fe/0 


3.29 


4.17 


2.79 


1.87 



Note. 



Baseline model corresponds to a coronal peak Alfven wave amplitude of ^ 60 km s ^, and 



the standard chromospheric model. ^ gives results with chromosphere de nsity modified to be co nsistent 
with photoionization-recombination equilibrium and the temperature in the lAvrett fc Loeserl (j2008[ ) model. 

give results with the ionization fract ion for each element m odified to be given by Saha equilibrium at 
the temperature and den s ity gi ven by Avrett fc Loeser ( 2008h . The final column shows the effect of the 
approximation of Laming ( 20091 ). where isotropic turbulence is assumed. 



Table 2: Open Magnetic Field Wave Spectra at 500,000 km Altitude 



ang. freq. 


vO 


vl v2 


0.010 


12.5 


12.5 125 


0.031 


150 


150 150 


0.062 


75 


75 75 


0.093 


50 


50 50 


0.124 


12.5 


125 12.5 



Note. 

Fig. 3 in lCranmer et al 



Frequencies are ang ular frequencies in rad s 
(l2007l) . 



Velocities are in km s "'^ ; vO is designed to match 
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Table 3: FIP Fractionations in Open Magnetic Field 







models 








observations 




ratio 


vO 


vl 




v2 


a 


b 




c 


He/0 


0.90 


0.85 





85 


0.60-0.58 


0.37-0.47 


0.45-0.55 


C/0 


1.13 


1.18 


1 


10 


1.50-1.41 


1.17-1.35 


0.9 


- 1.1 


N/0 


0.96 


0.94 





94 


1.19-0.9 


0.64-0.99 






Ne/0 


0.95 


0.92 





92 


0.48-0.47 


0.40-0.56 


0.3 


- 0.4 


Na/0 


1.97 


2.99 


2 


04 










Mg/0 


1.65 


2.21 


1 


67 


1.73-1.92 


0.98-1.60 


0.95 - 


2.45 


Al/0 


1.72 


2.37 


1 


75 










Si/0 


1.51 


1.89 


1 


53 


2.07-1.92 


1.20-2.09 


0.9 


- 1.8 


S/0 


1.26 


1.39 


1 


26 


1.53-1.56 


1.38-2.57 






Ar/0 


1.00 


0.99 





99 










K/0 


2.03 


3.14 


2 


12 










Ca/0 


1.88 


2.74 


1 


93 










Fe/0 


1.97 


2.97 


2 


04 


1.42-1.73 


1.04-1.69 


0.65 - 


1.35 


Ni/0 


1.85 


2.67 


1 


91 










Kr/0 


1.01 


1.01 


1 


01 










Rb/0 


1.97 


2.94 


2 


04 










W/0 


1.99 


2.99 


2 


07 











Note. — Prom left to right, model FIP fractionations correspond to the chromospheric model in Pig. 3. 

Observational ratios are taken from, (a) Zurbuchcn et al. (2002), (b) von Steiger et al. (2000), (c) Ko ct al. 
(2006), all given relative to O. Ranges quoted from von Steiger et al. (2000) include their uncertainties. 
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Table 4: FIP Fractionations in Closed Magnetic Field I 







models 








observations 








0.4 


0.5 


0.6 


a 


b 


c 


d 


e 


f 


ratio 


(km s ^ 


) 














He/0 


0.72 


0.61 


0.47 


0.68-0.60 


0.29-0.75 










C/0 


1.01 


1.03 


0.97 


1.36-1.41 


1.06-1.37 










N/0 


0.86 


0.81 


0.69 


0.72-1.32 


0.22-0.89 










Ne/0 


0.86 


0.71 


0.57 


0.58 


0.38-0.75 










Na/0 


2.42 


3.96 


6.74 






7.8+- 


1.8+2 






Mg/0 


1.76 


2.35 


3.25 


2.58-2.61 


1.08-2.36 


2.8!;;;5 


2.7±0.3 






Al/0 


1.95 


2.82 


4.12 








K «+3-3 

5-6_2.i 






Si/0 


1.70 


2.27 


3.02 


2.49-3.11 


1.36-3.24 


4.9 «;» 








S/0 


1.34 


1.57 


1.78 


1.62-1.92 


1.23-2.68 


2.2±0.2 


2.1±0.2 


1.7 ±0.3 




Ar/0 


0.96 


0.94 


0.86 










1.1 ±0.1 


1.12±0.15 


K/0 


2.72 


4.70 


8.51 






1 S+O-4 


4 7+7-0 

^- ' -2.8 


3.5 ±0.9 


6 


Ca/0 


2.38 


3.80 


6.27 






o i;+4.3 
"^■^-1.9 


2.7±0.25 




3.0-9.7 


Fe/0 


2.65 


4.44 


7.85 


2.28-2.90 


0.96-2.46 


7 0+1-4 

' -"-1.2 








Ni/0 


2.59 


4.10 


6.92 














Kr/0 


0.99 


1.00 


0.92 














Rb/0 


2.84 


4.96 


9.01 














W/0 


3.01 


5.35 


9.85 















Note. — Prom left to right, model FIP fractionations correspond to the chromospheric model in Pig. 3. 

The models for 0.4, 0.5 and 0.6 km refer to the amplitude with which the coronal wave is initiated in 
the low chromosphere. The corresponding coronal wave amplitudes are 55, 67, and 82 km s~^ respectively. 
Observational ratios are taken from, (a) Zurbuchen et al. (2002), given relative to O, (b) von Steiger et al. 
(2000), relative to O, (c) Bryans et al. (2009), given relative to the mean of O, Ne and Ar, (d) Giammanco 
et al. (2007, 2008), relative to H, and (e) Phillips et al. (2003), relative to H, (f) Sylwester et al. (2010ab) 
and McKenzie & Feldman (1994), relative to H, all given relative to the photospheric abundance of Grevesse 
& Sauval (1998). Ranges quoted from von Steiger et al. (2000) include their uncertainties. 
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Tablc 5: FIP Fractioiiatioiis in Closed Magnetic Field II 



ang. freq. (rad s ^) 


U.Uoo 


U.Ud 


U.UbO 


U.U/ 


U.U/ 


n HQ 
U.Uo 


U.Uoo 


u.uy 


U.iUo 


Vinit (km s~^) 


U.ZOO 


U.zo 


U.oo 


U.4o 


U.OO 


U.O 


U.4o 


U.OO 


r\ 1 no 


Vcor (km s"-^) 


A n 


4U 


A c; 
4o 






dU 


oU 


OO 
OO 


on 
zU 


ratio 




















He/U 


U.o4 


0.81 


(J. 74 


0.63 


0.60 


U.74 


0.86 


0.94 


n no 

0.98 


U/U 


1 AO 


1.38 


1 oo 

l.oz 


1.18 


1.26 


1 /in 
1.4U 


1.53 


1.59 


1.56 


iN/U 


u.yo 


0.92 


0.89 


0.83 


0.80 


0.88 


n n /I 

0.94 


n nT 

0.97 


n nn 

0.99 


iNe/ w 


u.yu 


u.oy 


U.o4 


U. / 


yj. 1 z 


u.oo 


n Qi 
u.yi 


u.yo 


u.yo 


Na/0 


4.74 


4.47 


4.61 


4.31 


4.47 


4.36 


4.69 


4.73 


4.60 


Mg/0 


3.29 


3.13 


3.11 


2.74 


2.89 


3.13 


3.49 


3.59 


3.51 


Al/0 


3.59 


3.43 


3.49 


3.15 


3.30 


3.42 


3.73 


3.81 


3.75 


Si/0 


2.70 


2.63 


2.69 


2.51 


2.69 


2.74 


2.88 


2.89 


2.84 


S/0 


1.79 


1.77 


1.79 


1.72 


1.83 


1.86 


1.92 


1.90 


1.87 


Ar/0 


0.98 


0.98 


0.97 


0.96 


0.94 


0.97 


0.99 


0.99 


1.00 


K/0 


5.08 


4.85 


5.13 


4.84 


4.92 


4.68 


4.95 


4.98 


4.94 


Ca/0 


4.31 


4.13 


4.32 


4.00 


4.12 


4.06 


4.34 


4.39 


4.36 


Fe/0 


4.78 


4.60 


4.87 


4.52 


4.57 


4.44 


4.73 


4.79 


4.79 


Ni/0 


4.20 


4.09 


4.35 


4.12 


4.23 


4.06 


4.24 


4.25 


4.27 


Kr/0 


1.00 


1.00 


1.00 


1.00 


0.99 


1.00 


1.00 


1.00 


1.00 


Rb/0 


4.73 


4.61 


4.98 


4.79 


4.80 


4.50 


4.67 


4.67 


4.72 


W/0 


4.84 


4.73 


5.15 


4.83 


4.83 


4.56 


4.76 


4.79 


4.89 



Note. — Entries in rows Vinit and Vcor refer to the amplitude with which the coronal wave is initiated in 
the low chromosphere, and it corresponding amplitude in the corona. 
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corona 




Fig. 1. — Cartoon illustrating the model. Alfven waves generated in the corona either reflect 
from each footpoint or precipitate down, depending on their frequency with respect to the 
loop resonance. Resonant waves reflect, nonresonant waves precipitate down. Upcoming 
waves in the chromosphere, deriving from the mode or parametric conversion of p-modes 
at the /3 = 1.2 layer, are generally reflected back downwards, as illustrated at footpoint 
B. In our models, we specify wave amplitudes at footpoint A, and integrate the non-WKB 
transport equation back to footpoint B, where FIP fractionations are evaluated. 
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Fig. 2. — Coronal section of loop, length 100,000 km, magnetic field 20 G, (half wavelength 
long for a 0.07 rad angular frequency Alfven wave) showing from top: Elsasser variables 
in km s~^ {8B j ^/^Kp solid lines, 8v dashed lines), with black lines for real parts and gray 
lines for imaginary parts. Middle; wave energy fluxes in ergs cm~^ s~^, the thin solid line 
shows the difference in energy fluxes divided by the magnetic fleld strength and should be 
a horizontal line if energy is properly conserved. Bottom, the ponderomotive acceleration 
in cm s~^. Positive acceleration means positive along the z axis, which is upwards pointing 
near z — ^ and downwards near z — 100, 000. 
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Fig. 3. — Same as figure 2 giving the first three panels for the left hand chromosphere 
"B", where waves leak down from the corona. The extra bottom right panel shows the 
FIP fractionations (in black) for the ratios Fe/0, Mg/0, S/Oand He/0. Chromospheric 
ionization fractions are also shown in the fourth panel (in gray, to be read on the left hand 
axis) in the same hnestyles as for the fractionation. The extra dash-triple dot gives the 
O ionization fraction, the long dash curve gives the H ionization fraction. Fe and Mg are 
essentially fully ionized throughout the fractionation region. 
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Fig. 4. — Same as figure 2 for an open field region. Five waves corresponding to the baseline 
model in Table 2 are considered, initiated at 500,000 km altitude and integrated back to the 
chromosphere. 
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Fig. 5. — Chromosphcric section of open field calculation of Figure 4. Fractionations are 
shown for Fe/0, Ar/0 and He/0. 
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Fig. 6. — Ponderomotive force (top panels) and FIP fractionations (bottom panels) for the 
chromosphere including five chromospheric waves and a coronal wave with angular frequency 
0.06 (left) and 0.07 rad s~^ (right) respectively. Fractionations for Fe/0, Mg/O, S/0, and 
He/0 are shown. 
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Fig. 7. — Ponderomotive force (top panels) and FIP fractionations (bottom panels) for the 
chromosphere including five chromospheric waves and a coronal wave with angular frequency 
0.085 (left) and 0.0105 rad s"^ (right) respectively. Fractionations for Fe/0, Mg/0, S/0, 
and He/0 are shown. The 0.085 coronal wave gives a more similar FIP effect for Fe/0 and 
Mg/0, as frequently observed. The He/0 depletion reduces as the coronal wave moves off 
resonance, as the "spike" in the ponderomotive acceleration decreases in prominence. 
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Fig. 8. — Illustration of FIP fractionation for the 0.07 rad wave, showing the corre- 
spondence with important features of the chromosphere. Top left gives the chromospheric 
density and temperature with height. Bottom left gives the ponderomotive acceleration as 
before, with solid and dashed curves showing the acceleration with and without energy loss 
to slow mode waves. The dotted curve shows the slow mode wave amphtude. The left panels 
show the "spike" in the ponderomotive acceleration at the altitude where the chromospheric 
density gradient is strongest. Bottom right is the same as before, with FIP fractionations 
for Fe/0, Mg/0, He/0, and S/0, together with ionization fractions. Top right shows the 
ionization fractions for C. S, Mg, and Fc in an expanded view. Thick curves correspond to 
the "baseline" charge state fractions used for FIP fractionation throughout this paper, thin 
curves give the results of the Saha approximation. This overestimates ionization fraction 
close to the top of the chromosphere, because photons resulting from radiative recombina- 
tions are not allowed to escape, but remain trapped to cause further photoionization. 



